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Abstract 

Singular spectrum analysis (SSA) is a nonparametric and adaptive 
spectral decomposition of a time series. The singular value decompo¬ 
sition of the trajectory matrix and the anti-diagonal averaging leads to 
a time-series decomposition. In this algorithm, a single free parame¬ 
ter, window length K, is involved which is the FIR filter length for the 
time series. There are no generally accepted criterion for the proper 
choice of the window length K. Moreover, the proper window length 
depends on the specific problem which we are interested in. Thus, it is 
important to monitor the spectral structure of the SSA decomposition 
and its window length dependence in detail for the practical applica¬ 
tion. In this paper, based on the filtering interpretation of SSA, it is 
shown that the decomposition of the power spectrum for the original 
time series is possible with the filters constructed from the eigenvec¬ 
tors of the lagged-covariance matrix. With this, we can obtain insights 
into the spectral structure of the SSA decomposition and it helps us 
for the proper choice of the window length in the practical application 
of SSA. 
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1 Introduction 


Singular Spectrum Analysis (SSA) is a nonparametric and adaptive decom¬ 
position of a time series. In the SSA, the original time series is separated 
into arbitrary number of additive components which are grouped and are 
interpreted as the slowly varying trend, oscillatory and noise components, 
respectively. It has been applied to many different fields from physics and 
earth sciences to financial and social sciences [2-7]. 

In the conventional treatment of SSA, the algorithm consists of four steps: 

• 1st step: a time series {xq,Xi,x 2 , ...,Xn-i} of length N is transformed 
into multidimensional series by constructing the L x K trajectory ma¬ 
trix X. Where K is the window length which is an arbitrary positive 
integer satisfying K < N/2 and L — N — K + 1. 

• 2nd step: The singular value decomposition of the trajectory matrix 

X = X 1 + X 2 + ... + X K . (1) 

• 3rd step: so-called grouping procedure where the above set of rank-one 
matrices X k is partitioned into several groups. 

• 4th step: reconstruction process where the decomposition of the time 
series is completed with the anti-diagonal averaging of the grouped 
matrix X^ = X n + X,; 2 + Xj 3 + ... + X lrn which transform X^ into 
Hankel form. 

The major objective of the SSA algorithm is to separate the original time 
series into interpretable set of time series. A single parameter, window length 
K, is involved in the algorithm and the results of decomposition depend on 
the choice of this parameter K. These points are throughly discussed from 
the separability point of view [3,4], It is argued that the relatively large 
window lengths are preferable to achieve a better (weak) separation, but 
for the series with a complex structure, too large window length leads to 
undesirable decomposition or mixture of the series components. Then, how 
to choose the proper window length K is an important point for the practical 
application of SSA but is not easy to formulate in a general way. The choice 
of the window lengths have been done so far by investigating the behavior 
of the eigentriples (singular values, and left and right singular vectors of the 
trajectory matrix) or the w-correlation matrix. 
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In this paper, based on the filtering interpretation of SSA [1,8,10,12], the 
decomposition of the power spectrum for the original time series is shown to 
be possible and this gives us detailed information about the spectral structure 
of the SSA decomposition. In the filtering interpretation, the SSA algorithm 
can be seen as the generation of normalized and complete set of adaptive 
filters and the operation of these filters to the original data. The SSA de¬ 
composition is thus interpreted as the soft and the optimal partitioning of 
the Fourier space through the adaptively constructed filters. By utilizing the 
completeness property of the eigenvectors of the lagged-covariance matrix, 
the filters constructed with these eigenvectors are used for the decomposition 
of the power spectrum of an original time series. The decomposition of the 
power spectrum shows the detailed structure of the SSA decomposition and, 
with this, we can obtain useful information about the spectral structure of 
the SSA decomposition and this helps us to choose the proper window length 
K in the practical application. 

In the next Sec. 2, we briefly summarize the conventional treatment of 
the SSA algorithm. We also describe the filtering interpretation of the SSA 
and it is adopted for the decomposition of the power spectrum of original 
time series. In Sec. 3, in order to see how the decomposition of the power 
spectrum works, we have examined two cases: the schematic time series with 
the superposition of the pure harmonic components and the daily currency 
exchange rates of EUR/USD. The summary and conclusions are given in Sec. 
4. 

2 Algorithm for the singular spectrum anal¬ 
ysis 

In this section, we briefly describe the SSA algorithm. In the filtering in¬ 
terpretation, the SSA algorithm consists of two steps. First, optimal and 
adaptive construction of the normalized and complete set of filters through 
the eigenvectors of the lagged-covariance matrix. Second step is the operation 
of these filters to the original time series. In the conventional treatment of 
the SSA, these two procedures are mixed into the singular value decomposi¬ 
tion and the anti-diagonal averaging procedure. Because of the completeness 
property of the filters generated in SSA, not only the time series but also the 
power spectrum can be decomposed into arbitrary number of additive com- 
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ponents. In order to show this, we assume the periodic boundary condition 
for the time series to avoid subtleties coming from the starting or the end 
points of a time series. 


2.1 Conventional SSA algorithm 

Let us consider a real-valued time series {x 0 , aq, x 2 , ..., aqv-i} of sufficient 
length. First, we define the L x K trajectory matrix X as 


Xij Xj_|_j_2 (i 1)2,..., T, j 1, 2,..., K') 


( 2 ) 


where window length K is the arbitrary positive integer satisfying K < N/2 
and L = N—K+l. Next step is the singular value decomposition of the above 
trajectory matrix X. To do this, the lagged-covariance matrix C = X T X is 
introduced and the eigenvalue equation is solved 

CV (fc) = A fc v (fc) . (3) 

The eigenvalues A*, are ordered as the decreasing order of magnitude (Ai, A 2 , •••, Xk) 
and the corresponding eigenvectors are denoted as v^. We define the L di¬ 
mensional vector as w ^ = Xv^ / y/Xk. Then the L x K rank-one matrices 
X k are expressed as X k = y/\j t 'w^'v^ T . Then the singular value decompo¬ 
sition of the trajectory matrix is expressed as 


X = y x k = y X^ k)v(t)T 

k= 1 k=1 


(4) 


The third step is the partitioning of the matrices X k into several groups and 
is called as the grouping procedure. The fourth step is the reconstruction of 
the time series. For the grouped matrix X^ = X n +X i2 +X l3 + ... + X im , we 
carry out the anti-diagonal average which transform into Hankel form. 
Then the time series are reconstructed giving the SSA decomposition. 


2.2 Filtering interpretation of SSA algorithm 

As are discussed in Refs. [8,10-12], above mentioned SSA algorithm can 
be viewed as the two-step point-symmetric filtering of the time series. The 
filtering interpretation is much more transparent in Fourier space, and the 
normalization and the completeness properties of the Liters play a crucial role 
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in SSA decomposition. To see this, we define the discrete Fourier transform 
as 


x a 


1 L ~ 1 

“Tr H exp(2vri aj / L) Xj , 

v T j=Q 


(5) 


and 

#2° = exp(2vri aj/L) 

3=0 


( 6 ) 


where the components of the vector vF) is denoted as 
In above Eq. 

,(*) „,( fc ) 


v 


(*) „,(*) 


1 


Vo 


"> u K 


y 


d]) , we consider the Fourier transform of the series of length L 
FI no m \^ l - i - g re pi ace( } by 


(uj ,i >2 ,v ( K , 0, 0,0) and then the sum J2j=o 10 ‘cpiatcu uy ^ - =0 
The forward filtering Iv can be rewritten in the Fourier space as 


Xr 


v^x = ■yF )’ 1 

U — a Js a u a 


Xr 


(7) 


The above relation holds for the periodic series. Since the eigenvectors are 

real valued, their Fourier transform satisfies = v^*. In above equation, 
the vectors v F) are complex in general, and shift and distort the time series. 
The non-causal FIR filtering Eq. (J7J) corresponds to the SVD of the trajectory 
matrix in the conventional treatment of SSA. The last step of SSA is the anti¬ 
diagonal averaging (Hankelization) for the rank-one matrix Xf. and this can 
be expressed as the causal FIR filtering 

■Xa -» V^Xa ■ ( 8 ) 


Then the two-step filtering can be combined as a single real filtering 


x a 




X n 


1 


(9) 


and does not distort the time series [1,9,10,12], It should be important to 
note that the orthonormal and the completeness properties of the eigenvector 
vF) lead to the following relations for the filters || 2 [10,11] 

El «?’I7£ = i 

Q!—0 

and 

Eief/A^i 

k=l 


( 11 ) 


( 10 ) 
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An arbitrary orthonormal iT-dimensional vector set {v^} leads to the fil¬ 
ters satisfying the relations Eqs. (TT01) and (ITT]) . In the SSA, the specific 
iT-dimensional vectors v are adopted which maximize the norm of the L- 
dimensional vector Tv preserving the orthonormality of the vectors v, which 
leads to the eigenvalue equation Eq. (J3]). For the eigenvalues and the eigen¬ 
vectors for the lagged-covariance matrix C, the following relation holds 

A fc = v^ T X T Xv^ = ^|h( fc) | 2 |x Q | 2 . (12) 

a 

The SSA can be interpreted as an algorithm to determine the orthonormal 
K -dimensional vectors v’D which maximize the quantity |h a | 2 |a: a | 2 pre¬ 
serving the orthonormality of the vectors v. In this sense, the the SSA can 
be understood as the principal component analysis applied to the trajectory 
matrix constructed from the lagged vectors. The orthogonality of the eigen¬ 
vectors are preserved as the phase content for v£\ and then, it is lost for the 
filters \v^\ 2 . Instead of loosing the orthogonality, the filters \v^\ 2 are rea l 
valued and lead no phase shift or distortion for the time series [1,9,10,12], 
The SSA decomposition of the time series can simply be expressed as 


x a = (l^l 2 x a + |v£ 2) | 2 x a + ... + |hi K) | 2 x a )/K . (13) 

We denote the inverse Fourier transform of the component \v^\ 2 x a as 
Then, the SSA decomposition of the time series is expressed as 

x n = x ^ + x^ 2) H-f x^ . (14) 

For later convenience, we denote the filters as 

A 1 = K'fl 2 . (15) 

and the Eqs. (TTOT) and (fTTT) can be rewritten as 

Ed‘Vr = i . (is) 

a 

and 

'£Fi k) /K = 1 . (17) 

k 


From the filtering viewpoint, the SSA algorithm can be seen as the optimal 
generation of the filter set satisfying the relations Eqs. (fTOj) and (ITT]) , and the 
time-series decomposition can be done as in Eq. (1131) . From the completeness 
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relation Eq. m, we can see that the decomposition of the power spectrum 
of the time series | 2 is possible as 


|£«| 


2 


1 

K 




k =1 


2 


(18) 


Through the above decomposition of power spectrum, we can obtain the 
insights into the spectral structure of SSA decomposition. On the other 
hand, squared Fourier transform of the SSA decomposition x^ in Eq. (fl3l) 
corresponds to [ \v^] 2 x a ] 2 = |f)W| 4 |x a | 2 and their sum does not reproduce 
the power spectrum of the time series. In the next section., we will show 
that the above decomposition can be applied to analyze the details of the 
SSA decomposition and it can be used to monitor the spectral structure of 
the decomposition which helps us for the proper choice of the window length 
parameter K. 


3 Example 

To illustrate how the spectral decomposition of the power spectrum described 
in the last section works, we consider the schematic time series and the daily 
currency exchange rate EUR/USD. 

3.1 Schematic model 

First, we adopt the simple time series which is the superposition of the har¬ 
monic components with additive noise, 

4 

x n = Y^ aism(27m/bi) + e n , (19) 

i =1 

where N = 300 and e n the Gaussian noise N(0, 1). We adopt the parameters: 
cq = 1.44, a 2 = 1.04, a 3 = 1.04, a 4 = 0.8; b 3 = 29.9, b 2 = 9.7, b 3 = 
8.9 and fe 4 = 4.9. The above time series Eq. (QUl) and its power spectrum 
are shown in Fig. 1. Corresponding to the harmonic components, there 
are four peaks A~D in the power spectrum. As seen, the second B and 
the third C peaks are close and have similar strengths. In such cases, these 
peaks are mixed with each other and are difficult to separate in SSA as is 
discussed in Ref. [3]. To see the dependence on the choice of the window 
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Figure 1: The upper figure is the time series in Eq. (TT9|) and the lower is its 
power spectrum. 

length, we adopt the values K = 40, 80 and 120. The larger 15 eigenvalues 
are shown in Fig. 2. It is known, for the case with larger window lengths, 
each harmonic component corresponds to the degenerate eigenvalue pair [3,4] 
and the corresponding eigenvectors are approximately the sinusoid with the 
relative phase shift 7 t/ 2 [3,13]. As seen, paired degeneracy of the eigenvalue 
becomes to be apparent for larger window lengths. For, K = 120, we can see 
clear four degenerate pairs. Since the filters are defined as the squared Fourier 
transform of the eigenvector of the lagged-covariance matrix as in Eq. (J3j), the 
Fourier transform v^ of the eigenvector for the degenerate pair differ only in 
phase and leads to approximately the same filters F^ k \ These filters extract 
the corresponding harmonic components from the time series. For all cases 
K = 40, 80 and 120, the peaks A and D are clearly filtered separately since 
these are isolated and the strengths are different from the other prominent 
spectral peaks. For the two peaks B and C, however, they are close with 
each other and the strengths are also similar. Though we can see the clear 
pairing of the eigenvalues for K = 120 in Fig. 2, it is not clear to what 



































Figure 2: The larger 15 eigenvalues of the lagged-covariance matrix for the 
window lengths K = 40, 80 and 120, respectively. We can see that the eigen¬ 
values are paired for larger K, corresponding to the harmonic components 
in the time series. 

extent these peaks are separated in the SSA decomposition. To see this, we 
have shown in Fig. 3, the decomposition of the power spectrum in Eq. (fl8jh 
Since the filter pairs (F^\ F^) and (F^ 7 \ F^) clearly extract the peaks A 
and D, respectively, these are not shown in the figures. The components 
^a^\^a\ 2 ~ Fjf)\x a \ 2 are shown for the window lengths K = 40, 80 and 120, 
respectively. As seen, the filtered strengths are fragmented and the peak D 
component are mixed for smaller window length K = 40. For K = 80, the 
component of the peak D almost disappears and the four filters F^ ~ F^ 
extract both B and C components. 

For larger window length K = 120, which is close to the half of the length 
of the time series, the Elters F^ and F^ extract mainly the component C 
with small admixture of B component, while the filters F^ and F^ extract 
the component B with small admixture of C component. The vector compo¬ 
nents for the eigenvectors and v (2 ^ are the sinusoid with 7 t/ 2 phase shift. 
The components of the paired eigenvectors (v®,v^) or (v^,v^) are the 
amplitude modulated harmonics with relative phase shift 7 t/ 2, representing 
the mixing of B and C components. These structure is easy to see from the 
decomposition of the power spectrum. Instead of this, we can calculate the 
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Figure 3: The spectral decompositions \x a \ 2 ~ F^\x a \ 2 are shown for 
the window lengths K = 40, 80 and 120, respectively. From left to right, 
the figures represent the spectrum strengths F^\x a \ 2 (k = 3 ~ 6) denoted 
as < k >. We take the same scale for all figures. 


squared Fourier transform of the decomposition x^ as \F^x a \ 2 but these 
quantities do not be summed up to reproduce the total power spectrum. The 
decomposition of the power spectrum has an advantage that the sum with 
respect to a corresponds to the eigenvalues as in Eq. (fl2l) and are summed 
up to reproduce the power spectrum as in Eq. (TT5T) . 

In order to see the details of the distribution of the spectral decomposi¬ 
tion, we have calculated the strength for individual spectral peaks as 

4 fc) = E F^\x a \ 2 i = A,B,C,D , (20) 

qG [neighbor of the peak i] 


and define the relative strengths as 


3 = 1 


K 


rU) 


( 21 ) 
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Figure 4: The relative strengths (k — 1 ~ 11) are shown for the spectral 
peaks i=A,B,C and D for the window lengths K = 40, 80 and 120, respec¬ 
tively. 


The results of the calculated relative strengths R \ k ^ are shown in Fig. 4. As 
is expected, the strengths coming from the paired filters dominate for the 
isolated peaks A and D for almost all values K considered here. For K = 40, 
the peak D component is extracted with the filters and F^ 6 \ while for 
larger window lengths, it is extracted with the filters F 1 ' 1 ' 1 and F^ s \ For 
the peaks B and C, the strength distribution depends on the choice of the 
window length K as seen in the Fig. 4. For smaller values of K , two peaks 
B and C are considerably mixed with each other. For larger window length 
K , the peaks B and C are almost separated with the paired filters but still 
remains the small admixture of the components B and C. Thus, the SSA 
decomposition (x^jX^) and (x^ 7 \ x^) are the pure harmonics with lower 
and higher frequencies, respectively, while the decompositions x ^ ~ x ^ are 
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the admixture of two harmonics. 

For the case that the spectral peaks are well separated, the smaller win¬ 
dow length is sufficient to separate the peaks of the spectrum and the results 
of the decomposition does not strongly dependent on the choice of the window 
length K. For the cases, where the spectral peaks are close or even overlap 
with each other, the Fourier peaks are mixed and also the SSA decomposi¬ 
tion is sensitive to the choice of the window length parameter [3,4]. Thus, for 
the time series with the complex spectral structure, we should be careful for 
the choice of the window length parameter. For the proper SSA decomposi¬ 
tion, it is helpful to monitor the spectral content of the decomposition and 
its dependence on the window length K by examining the power-spectrum 
decomposition as described here. 
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Figure 5: The daily currency exchange rate EUR/USD for the period 2013- 
2014. The number of data points is N = 489. The upper figure is the 
original daily rate and the first component of the SSA decomposition with 
the window length Jt=40. The lower graph is the residual time series. 


3.2 Daily currency exchange rate EUR/USD 

Next, we have applied the spectral decomposition for the analysis of the 
daily currency exchange rate Euro vs U.S. dollar for the period 2013-2014. 
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Figure 6 : The power spectrum of the residual components of the currency 
exchange rate EUR/USD shown as the lower graph in Fig. 5. 

The data consist of At=489 points. The trend component of the series is 
separated by subtracting the first component of the SSA decomposition with 
the window length K = 40, which is called as the sequential SSA [3]. The 
original exchange rate and the first component of the SSA decomposition are 
shown in the upper graph in Fig. 5. The residual component is also shown 
as the lower graph, and its power spectrum is shown in Fig. 6 . For the 
spectrum of the residual series, there are several peaks and are close with 
each other. Then, in order to see to what extent these harmonic components 
are separated, we have carried out the SSA decomposition with somewhat 
large window length parameters K = 140, 160 and 200. The eigenvalues 
Ai ~ A 13 are shown in Fig. 7. 

In this case, the eigenvalues (Ai, A2) and (A3, A4) are paired corresponding 
to the harmonic components in the series. The eigenvalues (A5, Ae) seems 
to be paired but are close with the values for the pair (A7, As) and these 
eigenvalues do not exhibit isolated pairing. There are five or six prominent 
peaks in the power spectrum and are distinguished by the attached indices 
A-F as in the upper graph of Fig. 8 . In the lower three graphs of Fig. 8 , we 
have shown the decomposition of the power spectrum for the case K = 140. 
Since the extracted components with pair of filters (F^\ F^), (F^ 3 \F^) 
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Figure 7: The larger 13 eigenvalues of the lagged-covariance matrix for the 
residual components with the window lengths K— 120, 160 and 200, respec¬ 
tively. The first and the second eigenvalues Ai and A 2 dominate. 

and (F^\ F^) are almost the same and we cannot distinguish them in the 
figure, we have shown the sum of these components (F^ +F^)\x a 1 2 . ui 3| + 
F^>)\x a \ 2 and (F^+F^)\x a \ 2 in these figures. The pair of filters (F (1 ', F^) 
corresponds to the largest peaks C, while the and (F^ 5 \ F^) 

pairs extract the peak D and E components, respectively. The broader peak 
B consists of the components A 7 and As, but in this case, there is moderate 
admixture of dominant C component. The peaks A and B are close with 
each other and the strengths are the mixture of Ag and A 10 components. In 
order to see the dependence on the choice of the window length K, we have 
calculated the relative strengths in Eq. (l 2 Tj) for the six peaks A-F and the 
results are shown in Fig. 9 for K = 120, 140 and 160, respectively. For 
the dominant harmonic components C, D and E, the results do not strongly 
depend on the choice of window length K. For the weak components A, B 
and F, the details of the filtering are considerably dependent on the choice 
of the window length K. For the separation of the harmonic components, 
window length around K = 140 seems to be appropriate. 

In the SSA, the filters F^ are optimally constructed from the time series. 


K = 120 
K = 140 
K = 160 
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Figure 8: The spectral decompositions F^\x a \ 2 ~ F^\x a \ 2 are shown for 
the window length K = 140. The individual components \x a \ 2 are ab¬ 
breviated as < k >. For the lower two graphs, a factor 6 is multiplied. 


We can say that the SSA algorithm is the soft partitioning of the Fourier 
spectrum with the normalized filters. Roughly speaking, the average widths 
of these filters in Fourier space gradually decreases with the increase of the 
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Figure 9: The relative strengths B^' } (k — 1 ~ 15) are shown for the peaks 
i —A ~ F with the window lengths K = 120, 140 and 160. 

window length K extracting the hner structure of the spectrum, but the 
structure of these filters are strongly dependent on the power spectrum of the 
time series itself. Generally, for dominant eigenvalues, corresponding filter 
has a simple structure of prominent peak and oscillating tails similar to 
the squared sine function. For non-dominant eigenvalues, however, the filters 
have complicated and broader function of a. Thus, the SSA decomposition 
is considerably dependent on the choice of window length and it is difficult 
to find the general criterion for the proper choice of the window length K. 
Because of these, the spectrum decomposition with the filters F ^ are useful 
to monitor the SSA decomposition in detail and adjust the window length in 
a specific application of the time series. 
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4 Summary and conclusions 

In this paper, based on the filtering interpretation of SSA, the decomposi¬ 
tion of the power spectrum for the original time series is shown to be possible 
and this helps us for obtaining the detailed insights into the SSA decompo¬ 
sition and it also helps us to choose the proper window length parameter K 
for a specific application. In the filtering interpretation, the SSA algorithm 
has been seen as the generation of normalized and complete set of adaptive 
filters and the SSA decomposition of the time series can be done through 
the operation of these filters to the original data. The SSA decomposition 
is the soft and the optimal partitioning of the Fourier space through the 
adaptively constructed real-valued filters. By using the completeness prop¬ 
erty of the eigenvectors of the lagged-covariance matrix, not only the time 
series itself but also the power spectrum can be decomposed into additive 
K components with the filters. The decomposition of the spectrum clearly 
shows the detailed structure of the SSA decomposition. Roughly speaking, 
the widths of the filters slowly decrease with the increase of window 
length K. This suggests that the finer structure of the spectrum of the time 
series can be separated with larger window length. But the structure of the 
filters are fairly complex and the the dependence on the window length is not 
so simple as to set any general criterion for the proper choice of the window 
length K. Besides eigentriplc of the SSA decomposition (eigenvalues, right 
and left singular values of the trajectory matrix), the decomposition of the 
power spectrum with the filters is useful for monitoring how the SSA algo¬ 
rithm works for the time-series decomposition and for the proper choice of 
the window length K. 
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